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The main leitmotiv of the present talk is analysis of different topological objects such as strings and 
domain walls in QCD. As it is well known, the standard model in general (and QCD in particular) 
does not support any kind of such objects due to some simple topological arguments. However, the 
situation cardinally changes when one considers QCD in large N c limit or in the limit of large baryon 
density when the unique light rj field emerges. I discuss rj- domain walls and rj global strings which 
occur in both systems: in high density QCD and in QCD at large N c . I also discuss the effects 
of quantum anomalies in the presence of topological objects in high-density QCD. The anomaly 
induced interactions lead to a number of interesting phenomena (such as the induced currents along 
the strings) which may have phenomenological consequences observable in neutron stars. 
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This talk is based on the recent works [1-4] devoted to the study of different topological objects in QCD. Domain 
walls and strings are common examples of topological defects which are present in various field theories. Domain walls 
are configurations of fields related to the nontrivial mapping ttq(M), while topologically stable strings are due to the 
mapping wi(M) where M is manifold of degenerate vacuum states. While it is generally believed that there are no 
domain walls or strings in the Standard Model due to the triviality of the corresponding mappings, such objects do 
. exist in QCD at large N c and in the regime of high baryon densities, where the symmetries of QCD are broken. The 
^ ' main reason why these two (apparently very different ) regimes lead to the similar topological objects is very simple: 
. in both cases the very light rj field emerges. 

This talk is organized as follows. Section 2 will explain how some very interesting topological structures appear in 
the regime of high baryon density. Section 3 is devoted to the analysis of the similar structures in QCD at large N c . 
Finally, Section 4 is devoted to the quantum anomalies in the presence of topological objects in high-density QCD. 
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*^ HIGH DENSITY QCD 

JL . The rich phase diagram of QCD at high baryon density has attracted considerable attention recently, see review 
paper [5] for references. A typical example is QCD with three light flavors. The ground state [the color-flavor-locked 
(CFL) phase] was determined to break baryon number, chiral and U(1)a symmetries. The CFL phase is characterized 
by the condensation of diquark Cooper pairs. These pairs are antisymmetric in spin (a, /?), flavor and color (a, b) 
indices: 
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X al = e t3k e abc e af3 {qi b q k p c )*, Y al = e ij V 6c e d/3 (^ b g| c )*. (1) 

The following object, £ = XY* in contrast to X and Y, is a gauge-invariant order parameter. Furthermore £ carries 
a nonzero U(1)a charge. Indeed, under the U(1)a rotations q — ► e 475 "/ 2 g, the X, Y fields transform as X — > e~ la X, 
Y -► e ta Y, and therefore E -> e" 2M £. 



Effective lagrangian for the rj field. 

The color-superconducting ground state, in which (£) ^ 0, breaks the U(1)a symmetry. The Goldstone mode 77' 
of this symmetry breaking is described by the phase cpA of S = |£|e~ i¥ '' 4 . Under the U(1)a rotation tpA transforms 
as ifA — * Pa + 2a. At low energies, the dynamics of the Goldstone mode tpA is described by an effective Lagrangian, 
which, to leading order in derivatives, must take the following form, 

L = f[(docpA) 2 -u 2 (d i <p A ) 2 ]. (2) 

It is well known that the U(1)a symmetry is not a true symmetry of the quantum theory, even when quarks are 
massless. The violation of the U(1)a symmetry is due to nonperturbative effects of instantons. Since at large 
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chemical potentials the instanton density is suppressed (see below), the rf boson still exists but acquires a finite mass. 
In other words, the anomaly adds a potential energy term V[ ns t(p) to the Lagrangian (2), 

L = f[(d <p A ) 2 - U 2 (d lVA f] - V inst (p A ) • (3) 

The curvature of VJ ns t around pa = determines the mass of the rj' '. Also, the standard arguments suggest that 
the potential Vi nst must depend on the single combination, pa — 0. Moreover, at large p,, Vi nst can be found from 
instanton calculations explicitly. The infrared problem that plagues these calculations in vacuum disappears at large 
p: large instantons are suppressed due to Debye screening. As a result, most instantons have small size p ~ p~ x and 
the dilute instanton gas approximation becomes reliable. One-instanton contribution, proportional to cos(pa — 8), 
dominates in Vinst- Therefore, 

Vinst {'Pa) = -ap 2 A 2 cos(p A - 9) , (4) 

where A is the Bardeen-Cooper-Schricffcr (BCS) gap, and a is a dimcnsionlcss function of p which can be explicitly 
calculated [1]. Here we only note that a vanishes in the limit p — > oo very fast, a <~ ( c j t CD ) b , b = 11/3A C — 2/3 AT/ . 
This is an important fact, since it implies that the mass of the rf boson, 

m=V^72^A, (5) 

becomes much smaller than the gap A at large p. In this case the effective theory (3, 4) is reliable, since meson modes 
much heavier than rf decouple from the dynamics we are interested in. 

The Lagrangian (3) with the potential (4) is just the sine-Gordon model, in which there exist domain- wall solutions 
to the classical equations of motion. The profile of the wall parallel to the xy plane is pa = 4 arctanexp(mz/w), so 
the wall interpolates between p = 0(z = — oo) and ipA = 2n(z = oo). The tension of the domain wall is 

a = 8V2aufpA. (6) 

One should note that in the limit a = the spontaneously broken U(1)a symmetry leads to there existence of the 
global strings [3] which however are always accompanied by the domain walls for any small but non -vanishing a/0. 



Decay of the domain wall. 

It is important to understand the mechanism of the decay of the wall. It has nothing to do with the decay of 
r/ 1 meson quanta, which are due to r)' coupling to photons, ungapped quarks, or the gluons of the unbroken SU(2) C 
subgroup. The domain wall is already a local minimum of the energy, and the decay of its excitations means only 
that the fluctuations around this minimum, corresponding to deformations of the wall, are damped. 

The domain wall is not stable because the same ground state is on both of its sides: p>A — and pa = 2tt are 
equivalent. The instability is due to higher energy meson modes integrated out and not present in the Lagrangian (3). 
One can visualize the effect of these modes by considering an effective potential which depends on the magnitude |E| 
as well as on the phase p of the order parameter S. This potential has the shape of a Mexican hat, slightly tilted by an 
angle proportional to a. The domain wall is a configuration that, as a function of the coordinate perpendicular to the 
wall, starts from the global minimum, goes along the valley, and returns to the starting point. One can continuously 
deform this configuration into a trivial constant one by pulling the looplike trajectory over the top of the hat. This 
deformation has to be done in a finite area of the wall first, thus creating a hole. If this hole exceeds the critical size, 
it will expand, destroying the wall. On the rim of the hole the magnitude of |E| vanishes. The field configuration 
around the rim is a vortex: on a closed path around the rim, pa changes by 2tt. The decay of the wall is a quantum 
tunneling process in which a hole bounded by a closed circular vortex line is nucleated. The semiclassical probability 
of this process is[l], 

where v is the tension of the vortex line in the limit of massless boson, m = 0. 

To find r we need to compute the string tension is, which to the logarithmical accuracy can be estimated as 
v = 2iru 2 f 2 ln(i?A) , where R is a long-distance cutoff (thickness of the wall <~ m _1 ), and A is the momentum 
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scale at which the effective Lagrangian description breaks down. We are helped by the fact that the vortex tension 
is dominated by the region outside the core, so the effective Lagrangian (2) is sufficient for computing v to the 
logarithmic accuracy. Therefore, the decay rate can be estimated to be 

r-cxpf-^^ln 3 ^). (8) 
\ 3 a (j, 2 A 2 sja) 

Since / <~ /j » A, and a decreases with increasing /x, the decay rate is exponentially suppressed at high \x. 



QCD AT LARGE N c . 

In this section we want to demonstrate that all important ingredients of the rf string and domain wall configurations 
discussed in the previous section, remain in place for the case of QCD at large number of colors. Indeed, the most 
profound property of QCD at large N c , which is there existence of the light pseudo Goldstone rj' field, remains 
the same. This is due to the fact that new small parameter 1/N C <C 1 enters into the system replacing the small 
parameter Aqcd/h <C 1 discussed previously. Indeed, the 77' mass in the chiral limit m 2 ~ 1/N C vanishes[6]. This 
relation replaces the 77' mass formula (5) of the previous section. Similar to the previous case, the periodic properties 
of the 9 parameter in QCD along with the specific dependence of the potential on the unique combination (9 — ipa), 
imply there existence of the 77' domain wall. As in the previous case, the domain wall configuration corresponds to 
the path which interpolates between (fA — and tpA = 2-7T describing the same physical vacuum states with identical 
energies (this corresponds to the spontaneous breaking of the discrete symmetry, (fA = — > fA = 27r) . The width of 
the wall determined by small rj' mass. Similar to the previous case, the spontaneous (rather than explicit) violation 
of the continuous U(1)a symmetry gives rise to there existence of the 7/ string which is attached to the 77' domain 
wall. 

However, there is a difference between the present system and the system discussed in the previous section. In the 
case of dense quark matter all calculations were under complete theoretical control for sufficiently large /i. In the 
present case of QCD at large N c while a general iV c counting represents a very robust result, some specific properties 
of the potential which are responsible for 77' mass are model dependent. For a specific model for the 77' potential 
motivated by the 9 dependence in the form ~ cos(9/N c ) the 77' domain wall solution can be explicitly constructed and 
the corresponding wall tension can be explicitly calculated [2], 

' V "-/>,, (l - cos ~ ^ (9) 



N f - ,7T 1 \ 27V C , 

Similar to the previous section, the domain walls decay due to the the hole nucleation. To describe the corresponding 
process we should look for a corresponding "instanton" which is a solution of Euclidean (imaginary time, t = it) field 
equations, approaching the unperturbed wall solution at r — > ±00. In this case the probability T of creating a hole 
with radius R can be estimated as follows T ~ e~ 5 ° where So is the classical "instanton" action. If the radius of the 
nucleating hole is much greater than the wall thickness, we can use the thin-string and thin-wall approximation. (The 
critical radius R c will be estimated later and this approximation justified). In this case, the action for the string and 
for the wall are proportional to the corresponding worldsheet areas 

So = AnR 2 v- ^TrR 3 a. (10) 

The first term is the energy cost of forming a string: v is the string tension and 4irR 2 is its worldsheet area. The 
second term is energy gain by the hole over the domain wall: a is the wall tension and 4/37ri? 3 is its worldsheet volume. 
The world sheet of a static wall lying in the x-y plane is the three-dimensional hyperplane z = 0. In the instanton 
solution, this hyperplane has a "hole" which is bounded by the closed worldsheet of the string. Minimizing (10) with 
respect to R we find the critical radius and the action, 

R c = — , S = . (11) 

a 3a 2 

The obtained expression for So of course coincides with the one used in the previous section (7). To complete the 
calculations we need to estimate the string tension, which (as in the previous case) is mainly determined by the region 
outside the string core. The effective largangian description is justified in this region, and therefore, to logarithmic 
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accuracy one gets, v ~ jNff 2 , logi? ~ N c \ogN c , where we have set the outer radius of the string R is the same order 
as the thickness of the domain wall R <~ mT 1 <~ \[W C - Therefore, we have the following N c dependences 

v~N c logN c , a~y/W c , (12a) 

Rc~ V^clogTV c , m-^y/Nl, (12b) 

5-0 ~ N 2 C log 3 N c , r ~ e"^ 2 lo s 3 ^ . (12c) 



To conclude: the rf domain walls are large classical configurations with very long life-time r ~ cxp(N 2 log 3 N c ). We 
also notice that, although both the critical radius for nucleation R c and the wall thickness m _1 increase, the ratio 
mR c ~ log N c 3> 1 and the semiclassical approximation used to derive the decay rate (12) becomes justified in the 
large N c limit. 

The structure which is just described is very generic, and we expect that one should be able to recover this structure 
in the string holographic description of the r/ at large N c [7]. One can also argue that the qualitative picture described 
above remains valid even for physically relevant case N c = 3. In this case the domain walls will have very large life 
time [2], such that they can be even produced and studied at RHIC(relativistic heavy ion collider) [8] . 

QUANTUM ANOMALIES IN DENSE MATTER 

It is well known that anomalies have important implications for low-energy physics: the electromagnetic decay of 
neutral pions is a textbook example. In this section I describe the roles of anomalies in QCD at finite density. I 
follow our recent work with Son[4]. We shall see that topological objects discussed previously play a crucial role in 
this study. The effects coming from these anomalies are strikingly unusual; they reveal intricate interactions between 
the topological objects such as vortices and domain walls, the Nambu-Goldstone (NG) bosons and gauge fields. 

Anomalies 

One method to derive anomalous terms in effective theories is to put the system in external background gauge fields 
and require that the effective theory reproduces the anomalies of the microscopic theory. We thus consider QCD 
in the background of two U(l) fields: the electromagnetic field and a fictitious (spurion) field which couples 
to the baryon current. Only at the end of the calculations we will put = /in M , = (1,0), corresponding to a 
finite baryon chemical potential. (The technique resembles the one used in Refs. [9, 10] in a similar context.) The 
Lagrangian describing the coupling of quarks with B v is 

J? B = h v {}d v + \B V ) (13) 
(the baryon charge of a quark is taken to be 1/3). The Lagrangian is invariant under U(1)b gauge transformations 

^ e W3^ ; B„ -> B„ + (14) 

Let us assume that the baryon number symmetry is spontaneously broken, and the condensate is formed. The low- 
energy dynamics then contains a NG boson ipB, which is the U(1)b phase of the condensate. The transformation 
property of ips can be found from Eq. (14), 

<Pb -> <pb + 2/3 (15) 

where the factor 2 is the baryon charge of the U(1)b breaking order parameter (assumed to have the baryon charge 
of a Cooper pair of two baryons). We also assume the there exist a neutral NG boson which comes from breaking of 
a chiral symmetry. This can be a ir°, i] or f/ boson. To keep our discussion general, we introduce the current which 
creates this boson, 
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This current generates a chiral transformation, tpi — ► exp^aQf J^ipi ■ The NG boson is characterized by a phase 
ipA, which transforms under this transformation as (fiA ~^ + <Zac* , where is a number characterizing the axial 
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charge of the condensate. The factor | was introduced in Eqs. (16) so that the U(1)a charge of the chiral condensate 
(<Zl<7r) is one. In color superconducting phases, the current will be normalized so that q A = 2. 

The effective low-energy description must respect the U(1)b gauge symmetry and be invariant under the gauge 
transformations (14) and (15). Therefore, in the effective Lagrangian the derivative d^fB should always appear in 
conjunction with to make a covariant derivative D^Lp B = d^Lp B — 2B^. Vice versa, each occurrence of has to 
accompanied by a — ^d /1 ip B . 

In the presence of the background electromagnetic and U(1)b fields, the conservation of (16) is violated by triangle 
anomalies: 

= - Y^2 (e 2 C A77 F^F^ - 2eC ABl B^F tllJ + C ABB B^B^) (17) 

where F^ u = d^A v - d v A^ and B^ = d^B^ - d v B^; F^ = \e llva pF af) , B^ v = ^e^ ua/3 B af3 . The coefficients C's in 
Eq. (17) are given by 

Ca 77 = Qi(Qi) 2 , Cab 7 = QjQi, C'abb = ^ Q\- (18) 

The effective theory has to reproduce the anomaly relation (17). In the effective theory d^jf^ can be found from the 
change of the action under the chiral transformations: SS = — Jd 4 x d^aj^ = Jd 4 xad^j^. From the condition of 
anomaly matching one can deduce that the effective Lagrangian contains the following terms [4], 

J^anom - ^—d^A \e 2 C A ^A u F^ - 2eC A B 7 (jm„ - \d v ip B ) F^ - ±C A bb e» vaf3 (^n u - \d v ^ B ) d a d^ B ] (19) 
57T q A L 

The term proportional to Ca 77 in Eq. (19) describes two-photon decays like ir° — + 2j and rj' — > 2j. Such processes 
occur already in the vacuum. We do not consider these terms here. The new terms are the ones proportional to 
Cab 7 and Cabb- At first sight, these new terms either vanish identically or are full derivatives and cannot have any 
physical effect. This is true when the NG fields ip A and ip B are small fluctuations from zero. In particular, ir° does 
not decay to a phonon and a photon. However, as (p A , B are periodic variables, the action can be nonzero if cither or 
both variables make a full 2ir rotation. This occurs in the presence of topological defects like vortices or domain walls 
discussed in section 2. 



Anomalies in CFL phase 

Our discussion so far has been rather general. We now specialize ourselves on specific color superconducting phase, 
CFL phase, and discuss the physical consequences in the presence of topological defects in this phase. In the CFL 
phase, the gauge-invariant order parameters are constructed from diquarks, see eq.(l), as £ = X^Y which is a 3 x 3 
matrix which breaks chiral and U(1)a symmetries. Obviously, the CFL phase has baryon vortices, see Ref. [3] for 
details. In fact, a rotating CFL core of a star has to be threaded by a vortex lattice. There are also U(1)a domain 
walls[l] and U(1)a strings attached to them(as discussed in Section 2). In what follows we consider only few interesting 
phenomena which by no means exhaust all possible effects(related to the anomalous lagrangian (19)). Our choice 
simply reflects the author's preferences and the page constraint requirements for the present Proceedings. 

Axial current on a superfluid vortex. — To start, let us consider a superfluid (baryon) vortex. Let the vortex be 
parallel to the z axis and located at x = y = 0. Around the vortex ip B changes by 2n 7 and in the center of the vortex 
it is ill-defined. Outside the core e^ va P d a dp<p B — 0, and the term in Jz? a non that contains ip B vanishes. However, at 
the vortex core, as usual, (d x d y — d y d x )ip B — 2n6 2 (x±). The action becomes 

S'anom = ~7T dt dz d z ip A (20) 



12tt _ 

where the integral is a linear integral along the vortex line. Now let us recall that the axial current is obtained, 
from Nocther's theorem, by differentiating the action with respect to d^A- One sees immediately that there is an 
axial current running on the superfluid vortex, with magnitude 

JA ^ vk' < 21 > 



This implies that the system induces a strong classical weak (rather than electromagnetic) field along the vortex. 
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Magnetization of axial domain walls. — Let us consider an axial domain wall in an external magnetic field. Such 
domain walls exist at very high densities where instanton effects are suppressed as discussed in section 2. When the 
baryon field B v is treated as a background, B v = (/i, 0), the following term is present in the anomaly Lagrangian: 

^AB 7 = ^g^B.V^ A (22) 

where we have set = 2 and B is the magnetic field. Consider now a U(1)a domain wall streched along the xy 
directions. On the wall tp& has a jump by 2ir. Now turn on a magnetic field perpendicular to the wall, i.e., along the 
z direction. Equation (22) implies that the energy is changed by a quantity proportional to BS, where S is the area 
of the domain wall. This means that the domain wall is magnetized, with a finite magnetic moment per unit area 
equal to eC AB7 t 1 1 '(4tt) . The magnetic moment is directed perpendicularly to the domain wall. For CFL phase the 
magnetic moment per unit area is e/i/(67r). 

Currents on axial vortices . — The same effect can be looked at from a different perspective. One rewrites Eq. (22) 
into the following form, 

-S^anom = ^ijk-^-i djdkf A (23) 

Since CijkdjdkfA ~ 2ttS 2 (x±) on the vortex core, the action can be written as a line integral along the vortex, 

S anom = e -^^jdi.A (24) 

which means that there is an electric current running along the core of the axial vortex, which is similar to the axial 
current on baryon vortices. The magnitude of the current is given by 

Naturally, the electromagnetic current running along a closed vortex loop generates a magnetic moment equal to \jS, 
where S is the area of the surface enclosed by the loop. This is exactly what we found in Eq. (22). A large vortex 
loop, therefore, has a magnetic moment that can be interpreted as created by the current running along the loop, or 
as the total magnetization of the domain wall stretched on the loop. The two pictures come from the same term in 
the anomaly Lagrangian. The currents flowing on vortices are reminiscent of Witten's superconducting strings [11]. 

Axial Current in the Magnetic Field Background B. — The effective lagrangian (22) leads to another interesting 
phenomenon [13]. Namely, in the quark dense matter in the presence of the magnetic field B the effective lagrangian 
(22) implies that there is the axial current running along B with the corresponding current density proportional to 
the magnetic field, j A — eC ^ il B. If the system were a type II superconductor, one could calculate the total current 
J A = J j A dS carried by a fundamental magnetic flux J BdS = -^2tt. The obtained current 

jA = CWyM (2g) 

is analogous to the one calculated previously on the supcrfluid vortex (21). However, the currents on the magnetic 
fluxes (26) could be much more important phcnomcnologically than the currents on the supcrfluid vortices (21) 
because the density of the magnetic fluxes in neutron stars could be many orders of magnitude larger than density of 
superfluid vortices. 

Angular momentum on the axial vortices. — An axial vortex loop also carries angular momentum. To see that, let 
us place a vortex loop inside a rotating quark matter. We thus take B v — \iv v where v v is the local velocity and 
diVj — djVi = 2eijkUJk- The ABB term in the anomalous Lagrangian leads to the following expression, 



^ABB = • V^A (27) 



CabbM 2 
8tt 2 

which implies that the domain wall is characterized by a constant angular momentum per unit area, equal to 
CabbM 2 /(47t). As in the case of the magnetic moment, this angular momentum can alternatively be interpreted 
as being carried by an energy (mass) current running along the vortex surrounding the wall, and the value of the 
energy current is 

;™ = ^g^ = |X> (28) 
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The last representation is particularly illuminating, since it shows that the same current (25) that creates the magnetic 
moment also creates angular momentum. We can visualize this current as carried by quarks with different flavors /, 
where each quark has electric charge Qf and energy equal to the quark chemical potential /x/3. 

Future Directions 

We have seen that the quantum anomalies, especially those involving the baryon or axial current, lead to new and 
extremely unusual effects in high-density matter. Some effects may have consequences for the physics of compact 
objects, which are to be explored. Of particular interest is the finding that baryon vortices (and magnetic fluxes) 
carry axial neutral current, which interact with neutrinos and may affect the cooling process of neutron stars or/and 
may influence the neutron star kicks. 

All results in this section have been derived from the anomaly terms of the effective Lagrangian and do not rely 
on details of the microscopic QCD Lagrangian. It would be very interesting to understand the microscopic origin of 
the vortex currents found here. The first step in this direction was undertaken quite recently[12] where some explicit 
calculations on the microscopical level have been carried out in a simple model. In particular, it was shown that in 
some cases the results derived from the anomaly terms are truly topological and do not get any corrections even for 
finite temperature T, while in other cases the corrections could be of order one. 

All credits for this project should go to my close collaborators M. M. Forbes, M. Metlitski, E. V. Shuryak, D. T. Son 
and M. A. Stephanov. They are not, however, responsible for any typos, errors and mistakes unavoidably present in 
this text. This work was supported, in part, by the Natural Sciences and Engineering Research Council of Canada. 
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